Primary tunnel junction thermometry 
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We describe the concept and experimental demonstration of primary thermometry based on a four 
probe measurement of a single tunnel junction embedded within four arrays of junctions. We show 
that in this configuration random sample specific and environment-related errors can be avoided. 
This method relates temperature directly to Boltzmann constant, which will form the basis of the 
definition of temperature and realization of official temperature scales in the future. 

PACS numbers: 



Temperature is a relatively poorly known quantity in 
modern metrology. It is well recognized that the way the 
international temperature scale is currently realized, in 
particular towards low temperatures, needs to be seri- 
ously reconsidered. It is currently based largely on arte- 
facts which should be replaced by methods relating to 
thermodynamic temperature via Boltzmann constant ks 
[l|, 0]. Methods based on solid state tunnel junctions, 
Coulomb blockade thermometry (CBT) [j, |4J and shot 
noise thermometry (SNT) [a, IfJ , have both shown great 
promise as fc^ -based thermometers for metrology. How- 
ever, both of them fall short up to now, when it comes 
to sufficient absolute accuracy. In case of SNT, the lim- 
itations are mainly of practical nature, and can possibly 
be overcome by a careful design of the sensor and the 
measurement set-up. For CBT, an uncontrolled error 
source is of more fundamental concern: CBT involves a 
measurement of a series connection of nominally identical 
junctions. The inevitable spread in junction parameters 
leads, however, to an error, which can usually be made 
small, but which limits the accuracy in particular when 
the average junction size is small @, @|. In this letter 
we introduce and demonstrate a method, single-junction 
thermometry, SJT, which combines the advantages of ba- 
sic CBT thermometry, but which avoids the parameter 
dispersion induced errors altogether. We show theoret- 
ically that the errors can then be efficiently suppressed, 
and demonstrate the operation in experiment. 

In Coulomb blockade thermometry an array of tun- 
nel junctions shows a drop in its differential conduc- 
tance around zero bias voltage, because of the influence 
of single-electron charging effects. The ideal operation 
regime of a CB thermometer is determined by the ratio of 
the single-electron charging energy, Eq — e 2 /2C, where 
C is the (average) junction capacitance, and the ther- 
mal energy fc^T at temperature T such that Ec <C HbT. 
The measured conductance peak [see Fig. [Ha)] has two 
important characteristics, its full voltage width at half 
minimum, Vi/2, and its normalized (by asymptotic con- 
ductance at large voltages, Gt) depth AG/Gt- Here the 
first one is given by b.AAksT/e per junction, and serves 
as the primary thermometer, provided the junctions in 
the sensor are mutually identical. The latter one is in- 
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FIG. 1: The single junction thermometer (SJT). (a) A typical 
conductance curve of a thermometer, (b) Schematics of the 
SJT. The white bounded areas are conductors, and the grey 
interconnects are tunnel junctions. 



versely proportional to T. SJT thermometry is based on 
the same principle as CBT but there the objective is to 
measure the conductance of a single tunnel junction, em- 
bedded in a four probe configuration through lines con- 
sisting of arrays of tunnel junctions, see Fig. [lib). In this 
topology, the advantageous protection from the influence 
of electromagnetic environment is achieved. At the same 
time, this configuration abolishes any requirement of a 
uniform structure, because only one junction is probed, 
and the rest of the junctions, indeed not necessarily iden- 
tical, act as an environment for this one. 

We separate the theoretical analysis into two parts. 
First we consider the case where the influence of the en- 
vironment beyond the junction array can be neglected. 
We show that the measurement is perfect in this case even 
for a non-uniform structure. Then we use this result as 
the seed for analyzing the influence of the dissipative en- 
vironment on the performance of the thermometer, and 
show that with sufficiently long junction arrays the ac- 
curacy can be maintained at the desired level. 

The tunnelling rate through a junction i in forward 
(+) or backward (— ) direction with normal conductors 
in thermal equilibrium is given by Q 
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where Rt,i is the junction resistance, j(x) — x/(l — 
e -x/k B T^ anc i £p± j g change in electrostatic energy 
in tunnelling. We may separate this energy change as 
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6F^ = ±eVi + SE^ lil where Vi is the (mean) voltage 
drop across junction i and SE^ . is the internal energy 
change associated with charging the capacitors of the ar- 
ray. In the present analysis we limit ourselves to the low- 
est order result in Ec/kgT, which is what yields the ba- 
sic results in thermometry. In this spirit, we then expand 
Y ,i{8F±) ~ r 0l< (±cV i )+r / o ii (±eV i )^, i - Analytic cor- 
rections for lower temperatures can be obtained readily 
by expanding up to higher orders, but they will not be 
considered here. The current E through junction i can be 
obtained as E = e £ { „ } a({n})[T 4SF+) - r ,i(5Fr)]. 
Here, cr({n}) is the occupation probability of the charge 
configuration {n} on the islands within the array. With 
these premises, and by using identities y(x) — 7(— x) = x 
and 7'(x) + ^'(—x) = 1, we obtain 



eR T ,i ^ 

-^ch,J- (2) 

The internal charging energy for each charge configura- 
tion {n} is given by E ch = ^ J2{n}( < ^~ 1 )hj n i n h wnere 
C _1 is the inverse capacitance matrix of the junction ar- 
ray. We have neglected the offset charges on the islands 
since in the high temperature regime the charge distri- 
bution is quasi-continuous, (Snf.) ^> 1 [7]. Let us denote 
the islands surrounding the junction i by L and R. Then 
for the relevant processes only ul changes into til =F 1 
and tir to ur ± 1 , whereas all the other charge numbers 
remain constant. There are two junction connections to 
islands L and R in SJT. Evaluating 8E~L . as the differ- 
ence of -E c h for the charge configurations before and after 
the tunnelling event, and making use of the properties 
E{n} a ({ n }) = 1 and J2{n} n k cr({n}) = for all k be- 
cause of the symmetry of E c h, we obtain the normalized 
conductance of junction i, (G/Gr)i = Rr,i§w- as 



(3) 



where S t = e^Cr 1 )^ + (C^rr - 2(C~ 1 ) M ], g(x) = 
e x [e x (x -2) + x + 2}/(e x - l) 3 and v % = eV t /k B T. The 
result of Eq. ^ is in fact the basis of the standard CBT 
formula in linear arrays of junctions. Equation tells 
that conductance of a single junction in SJT is accurate 
as a thermometer in any array of junctions: the magni- 
tude of conductance suppression depends on the distri- 
bution of junction sizes via the capacitance matrix (Si), 
but the temperature can be determined unambiguously 
from, e.g., the half width of g(vi), if Vi can be measured. 
Error-free measurement of Vi is indeed possible in the 
configuration of Fig. [ljb). This happens since the volt- 
age measurement via two arrays is typically performed 
using an amplifier with very large input impedance (in 
any case much larger than the resistance of the junction 
arrays). Then, essentially no current flows through these 



two arrays, and there is no voltage drop across them. 
Thus Vi is indeed the voltage seen by the amplifier. 

The analysis above applies to the case where the con- 
nection to the bias sources and signal amplifiers has zero 
impedance. In practice this is not the case, and the envi- 
ronment impedance introduces errors to Coulomb block- 
ade thermometry which can be suppressed by using long 
arrays of junctions Q. Similarly, one can realize the sin- 
gle junction measurement which avoids the errors by the 
embedding arrays. To analyze the remaining errors in 
this case quantitatively, we may write the tunnelling rates 
if instead of of Eq. Q) as 



r± = 
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1 (E')P(SF l ± - E')dE'. 



(4) 



Here, P(E) originates from the environment theory of 
single-electron tunnelling [13] , and it yields the probabil- 
ity (density) of electron to exchange energy E when it 
tunnels: positive (negative) E refers to energy emission 
(absorption) by electron. In the first analysis above we 
thus assumed P(E) = 8(E), i.e., the Dirac delta func- 
tion. The same steps as in the ideal dissipationless envi- 
ronment above lead now to 
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Here, h{x) = e x (e x - x - l)/(e x - if. For P(E) = 8(E), 
Eq. ([5]) reduces naturally to ([3]). Equation §5§ yields an 
easy way to evaluate the influence of environment even 
in complex circuits. 

In general P(E) is obtained from the phase-phase cor- 
relation function J(t) with P(E) = dt exp[J(i) + 
j-Et] fl2| . Here we assume that the array is uniform and 
embedded in a resistive environment with resistance R. 
We are interested in conductance of one junction within 
an array. J(t) can be written as 



J(*)=7r^{[cot(B)(l- e -H)-M 
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Here r = t/(R cq C cq ), B = h/(2k B TR cq C l 
R/N 2 and C cq = NC. For a linear array of N junctions 
we have 5i = ^pfr. The results for the corresponding 
SJT are identical to these upon replacing N by N' + 1 
in the expressions above. Here N 1 is the number of junc- 
tions in each of the four surrounding arrays. Above we 
have assumed that the stray capacitances are small. 

Figure [2] shows numerical results based on the anal- 
ysis above. Figure HJa) shows the width as normalized 
to the ideal width in delta-function environment for a 
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FIG. 2: The influence of environment on SJT thermometry, 
(a) R dependence of the half width of the conductance dips 
for N' = 1, 2, 3, 7, 15 (N = 2, 3, 4, 8, 16 in linear arrays) from 
top to bottom, (b) The normalized depth of the conductance 
drop, with the same parameters as in (a). At low values of 
R, N grows from bottom to up. (c), (d) The corresponding 
quantities on the logarithmic resistance scale for a single bare 
junction (N' = 0, N = 1). In these plots E c /k B T = 0.1. 



junction in arrays with varying length and as a function 
of R. It is clear that one needs to protect the junction 
by a long array, if accurate measurement of temperature 
using V1/2 is to be obtained. For N — 2, an error of 
about 10% in the range 100 < R < 500 fi is expected. 
The impedance of the environment at high frequencies is 
approximately Z env = \/e/fj,, determined by the permit- 
tivity e and permeability fi of the medium. For vacuum 
its value is ~ 377 f2, and for a circuit on silicon it is 
a few times smaller. Such increase of V1/2 due to en- 
vironment in short arrays is supported quantitatively by 
experiments, see, e.g., Ref. Q- The depth of the conduc- 
tance dip is shown in Fig. [2jb) : it depends strongly on R 
only in short arrays. We note that: (i) Since the environ- 
ment is never known precisely in the experiment, there 
is almost no way to correct theoretically for such errors: 
the only working strategy then is to suppress these errors 
precisely by embedding the measured junction in a long 
array, (ii) The effect of error suppression is essentially 
proportional to iV~ 2 (if k B TRC/h < N). Therefore, an 
array with N ~ 50 is in principle sufficient for measure- 
ments with 10 -4 absolute accuracy, (hi) Embedding a 
junction in a very resistive environment [1 01 ] instead of a 
junction array is not the best strategy in thermometry ei- 
ther, which is indicated by the very slowly decaying tails 
of the error at large values of R. To fully appreciate this 
point, we show in Fig. [21(c) and (d) the width and depth, 
respectively, of a single junction peak in a purely resis- 
tive environment. Although the width at large values of 
R slowly approaches unity, experimentally it is hard to 
fabricate resistive environments with R 3> 10 kfL 




FIG. 3: Samples and comparison of the data on SJT and 
a bare single junction, (a) Electron micrographs of the ref- 
erence structure with one junction connected to four leads 
(top), and the SJT structure with N' = 20 junctions in the 
leads (centre). Zoom of the central junction and of a section 
of a junction array are shown at the bottom, (b) Measure- 
ment of the conductance of the bare single junction and the 
SJT in Sample B at T ~ 0.3 K. The deeper and narrower 
drop in conductance corresponds to the SJT. The calculated 
conductance curves for the two samples based on the model 
described are shown by the solid lines assuming R = 80 Q,. 



This concludes our proof that, theoretically, the in- 
fluence of uncontrolled error sources, the inhomogeneity 
of the junction array and the noise of the environment, 
can be efficiently suppressed in a SJT. We discuss next 
the proof-of-the-concept experiments. Samples [see Fig. 
[3f a)] were fabricated by electron beam patterning and 
shadow angle evaporation with an oxidation step between 
the two electrode layers. Both the bottom and the top 
electrodes are of aluminium, they are 40 nm and 45 nm 
thick, respectively. The bottom electrode was thermally 
oxidized at 100 mbar for 10 min before deposition of the 
top electrode at an oblique angle. Two samples (A and B) 
with two types of structures have been measured in this 
work. The single tunnel junction was connected either 
directly to the external leads or it was embedded within 
four arrays of N' = 20 junctions, respectively. The two 
types of structures were fabricated on the same chip in 
the same vacuum cycle. Nominally, the central junctions 
are identical in the two cases, and all the junctions are 
0.6 fim 2 of area, yielding a junction resistance of ~ 6 kf2 
(Sample A) and ~ 4 kfi (Sample B). 

The samples were measured in a dilution refrigerator 
with a 40 mK base temperature. However, these struc- 
tures were not suitable for very low temperature mea- 
surements: we observe strong self-heating due to weak 
electron-phonon coupling in the present geometry near 
base temperature [13j . Therefore we present here data at 
temperatures at and above 150 mK. Conductance mea- 
surements in the SJT configuration yield a deeper and 
narrower peak than for the unprotected single junction 
structure in agreement with the calculated results of Fig. 
H This is shown in Fig. EDJb) for Sample B. Both the 
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FIG. 4: Measurements on a SJT (Sample A). Conductance 
dips at four temperatures are shown together with fits accord- 
ing to the model presented, including the self-heating correc- 
tion The inset shows a comparison of the temperature 
deduced from SJT (vertical axis) against that obtained by an 
ordinary CBT measurement across one of the embedding ar- 
rays (horizontal axis). The solid line has unit slope, and the 
error bars indicate the confidence interval of the fits. 



judge is the influence of the size of the islands between 
the junctions, I: as long as it is smaller than the dis- 
tance to "horizon", i -C Hc/kpiT, the above lump ele- 
ment analysis is valid 3, 17 1. Here, c = (^e) _1 is the 
signal propagation speed. The condition gets critical at 
particularly high temperatures, and extra care to place 
the array close to the junction has to be taken then. We 
want to add that the presented thermometry is not nec- 
essarily limited to the standard planar tunnel junction 
design, but may be applicable, e.g., in scanning probe or 
break junction geometries, since the junction parameters 
of the surrounding arrays need not be the same as those 
of the central junction. This might lead to the possibility 
of using tunable tunnel junctions in thermometry. 

Summarizing, we propose an absolute single tunnel 
junction thermometer. We have demonstrated the con- 
cept in preliminary experiments. Our method may turn 
out to be valuable in future realization of the interna- 
tional temperature scale based on Boltzmann constant. 

We thank Mikes (Centre for Metrology and Accredita- 
tion) and NanoSciERA project NanoFridge for financial 
support, Dmitri Averin for a useful discussion, and Mikko 
Mottonen for help in the analysis. 



SJT (N' = 20) and the bare junction (N' = 0) data fol- 
low the environment calculation assuming R = 80 fl, a 
value which is consistent with the discussion above. The 
arrays can thus indeed be employed to efficiently protect 
the junctions against environment fluctuations. 

Figure |4] shows four-probe conductance measurements 
of the SJT sample at a few temperatures together with 
fits according to the model presented. We included in the 
fits the influence of self-hcating in the manner presented 
in Ref. (HI, and this yielded perfect match to the peaks, 
see the lines on top of the data. The self-heating has, 
however, only a small influence on the curves at the tem- 
peratures shown. The extracted temperatures from such 
fits are shown in the inset of Fig. [4] against the reading 
of the CBT "reference thermometer", which was one of 
the N' — 20 junction arrays in the same sample. The 
dominating discrepancies between the two thermometers 
are due to finite errors in the fitting procedure in each 
case. The agreement between the two is good over the 
whole temperature range in Fig. 21 

In the theoretical analysis we focused on the high tem- 
perature and high junction resistance limit, and did not 
discuss errors due to, e.g., enhanced Coulomb effects at 
lower temperatures I8J and strong tunnelling in low resis- 
tance junctions [3, Ha). Yet these errors can be treated 
similarly to what has been done in standard Coulomb 
blockade thermometry, and their influence can be esti- 
mated and kept at a tolerable level by proper choice of 
junction sizes for each temperature and with suitable tun- 
nel barrier parameters. One more (controllable) error to 
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